arXiv:astro-ph/0211258vl 12 Nov 2002 


astro-ph/0211258 
CERN-TH/2002-207 
UMN-TH-2115/02 
TPI-MINN-02/44 


Updated Nucleosynthesis Constraints on Unstable Relic Particles 

Richard H. Cyburt^, John Ellis^, Brian D. Fields^ and Keith A. Olive^ 

^Department of Physics, University of Illinois, Urhana, IL 61801 
‘^Theory Division, CERN, CH 1211 Geneva 23, Switzerland 
^Department of Astronomy, University of Illinois, Urhana, IL 61801 
^Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

Abstract 

We revisit the upper limits on the abundance of unstable massive relic particles provided 
by the success of Big-Bang Nucleosynthesis calculations. We use the cosmic microwave 
background data to constrain the baryon-to-photon ratio, and incorporate an extensively 
updated compilation of cross sections into a new calculation of the network of reactions 
induced by electromagnetic showers that create and destroy the light elements deuterium, 
^He, ^He, ®Li and ^Li. We derive analytic approximations that complement and check the full 
numerical calculations. Considerations of the abundances of ^He and ®Li exclude exceptional 
regions of parameter space that would otherwise have been permitted by deuterium alone. 
We illustrate our results by applying them to massive gravitinos. If they weigh ~ 100 GeV, 
their primordial abundance should have been below about 10“^^ of the total entropy. This 
would imply an upper limit on the reheating temperature of a few times 10^ GeV, which 
could be a potential difficulty for some models of inflation. We discuss possible ways of 
evading this problem. 
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1 Introduction 


Some of the most stringent constraints on unstable massive particles are provided by their 
effects on the abundances of the light nuclei produced primordially in the early Universe via 
Big-Bang Nucleosynthesis (BBN) [1]-[14]. Astrophysical determinations of the abundances 
of D, ^He and ^Li agree with those found in homogeneous BBN calculations, for a suitable 
range of the baryon-to-photon ratio rj [15, 16]. The decay products of massive unstable 
particles such as gravitinos would have produced electromagnetic and/or hadronic showers 
in the early Universe, which could have either destroyed or created these nuclei, perturbing 
this concordance. Maintaining the concordance provides important upper limits on the 
abundances of such massive unstable particles. The power of this argument has recently 
been increased by observations of the power spectrum of fluctuations in the cosmic microwave 
background (CMB) [17]. These now provide an independent determination of rj that is in 
rather good agreement with the value suggested by BBN calculations [18]-[22], reducing one 
of the principal uncertainties in the previous BBN limits on massive unstable particles. 

This development has triggered us to re-evaluate these BBN constraints. We do so via a 
new calculation of the network of nuclear reactions induced by electromagnetic showers that 
create and destroy the light elements deuterium (D), ^He, ^He, ®Li and ^Li. We perform a 
coupled-channel analysis of the light-element abundances, which enables us to analyze the 
possible existence of isolated ‘islands’ of parameter space that are not minor perturbations 
of standard homogeneous BBN calculations. In carrying out this program, we also take the 
opportunity to improve previous theoretical treatments of some reactions in essential ways. 
We also derive analytic approximations to our results, which serve as a check of the numerics, 
and offer additional insight into the essential physics. 

Considerations of the D abundance alone would have led to the apparent existence of 
hnely-tuned ‘tails’ of parameter space that extend the region allowed by standard BBN 
calculations to rather larger abundances of unstable relic particles. In the past, these ‘tails’ 
were argued to be excluded on the basis of the combined D -|- ^He abundance. It is now 
widely considered that the complicated chemical and stellar history of ^He renders such 
an argument unsafe [23]. However, we show here that these tails are excluded robustly 
by the astrophysical abundance of ^He, and also conflict with the measured ®Li abundance. 
These other abundances also exclude a disconnected ‘channel’ of parameter space that would 
have been allowed by the D abundance alone. Overall, the most stringent upper limit on 
the possible abundance of an unstable massive relic particle X with lifetime tx ~ 10^ s is 
provided by the ®Li abundance. For tx = 10® s, we hud 

„o 

Mx^ < 5.0x10-^^ GeV, (1) 

with the upper limit from the ^He abundance being about two orders of magnitude less 
stringent. 

As an illustrative application of this new analysis, we reconsider the allowed abundance 
of unstable gravitinos G with masses up to about 10 TeV, which are expected to have 
lifetimes ^ 10^ s. Using standard calculations of thermal gravitino production in the early 
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Universe [24] - [32] , [8], the constraint (1) suggests a stringent upper limit on the reheating 
temperature Tr of the Universe following inflation: 

Tr < 10^^ GeV for ~ 100 GeV. (2) 

Such a stringent upper limit can be problematic for standard inflationary models [33], some 
of which predict higher Tr ~ 10^^ GeV. A low reheating temperature might also be problem¬ 
atic for some models of baryogenesis in which thermal production of baryon-number-violating 
particles is necessary, though not for others. The simplest leptogenesis scenarios [34] require 
the production of a right-handed neutrino state whose decays violate lepton number. How¬ 
ever, even scenarios which require the thermal production of right-handed neutrinos can in 
fact accommodate very low reheat temperatures [35, 36], and non-thermal production of 
particles with masses less than the inflaton mass is possible for very low reheating tempera¬ 
tures [37]. Leptogenesis scenarios [38] of this type as well as those involving preheating [39] 
have been proposed 

Reheat temperatures higher than (2) would be allowed if Tq ^ 10^ s, as might occur 
if Mq > 3 TeV and/or the gravitino has many available decay modes. Alternatively, one 
may consider the possibility of a very light gravitino that would be stable or metastable. 
However, in this case the upper limit (1) would apply to the next-to-lightest supersymmetric 
particle (NLSP). In conventional MSSM scenarios with a very light gravitino, this might be 
the lightest neutralino y. However, the relic y abundance may also be calculated, and is 
likely to conflict with the upper limit (1). 

This brief discussion serves to emphasize the importance of evaluating the cosmological 
upper limit on the possible primordial abundance of unstable relic particles, to which the 
bulk of this paper is devoted. In Section 2, we discuss the most relevant photodissociation 
and photoproduction processes and describe their implementation in a code to calculate the 
network of reactions creating and destroying light elements in the early Universe. A brief 
discussion of the observational constraints used is in Section 3. Analytical and numerical 
results for light-element abundances are presented in Section 4. Our main constraints on 
unstable particles are described in Section 5, where we compare and combine the upper 
limits obtained by considering different light nuclei. Finally, in Section 6 we discuss in more 
detail the implications of our results for cosmological gravitinos. Section 7 contains some 
comments on particles with longer (and shorter) lifetimes. In Appendix A, we compile and 
discuss the cross sections used in this analysis. 

2 Photon Injection and Abundance Evolution 

As an example of the constraints imposed on the abundance of a heavy metastable particle 
by observations of light elements, we consider the radiative decay of a massive particle (such 
as a gravitino) with lifetime > 10^ s. The energetic decay photon initiates an electromagnetic 
shower, which in turn initiates a network of nuclear interactions. The decays of some unstable 
heavy particles also initiate hadronic showers, which would provide an additional set of 

^We also recall that Affleck-Dine [40] leptogenesis [36, 41] can be accomplished at low Tr. 
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constraints. These are typically important for shorter lifetimes, r ^ 10^ s. However, we 
restrict our attention in this paper to the network of reactions induced by electromagnetic 
showers, and comment on very short (and very long) lifetimes in Section 7. 


2.1 The Initial Degraded Photon Spectrum 


At the epoch of interest to us here, tx ^ 10^ sec, the massive gravitino^ is very non- 
relativistic, and can be treated as if at rest with respect to the background. We assume 
that a gravitino of mass Mx decays into a photon ( 7 ) and a neutralino (y), each with their 
respective energies. 


Ery = 


- ml 


2M 


X 


and Eyy = 


+ m'l 


2M 


X 


( 3 ) 


In the limit Mx 3> the energies become almost equal: E^ ~ E^^ \Mx. 

The primary photon with injection energy Eq interacts with the background plasma and 
creates an electromagnetic cascade. The most rapid interactions in this cascade are pair 
production yybg —> e+e“ off of background photons, and inverse compton scattering. These 
processes rapidly redistribute the injected energy, and the nonthermal photon spectrum 
rapidly reaches a quasi-static equilibrium as discussed in [7, 9, 10]. The “zeroth generation” 
quasi-equilibrium photon energy spectrum is 

f 

P-,{E-t) = < A'o (|^)^° for Ex < E., < Ec 
{ 0 otherwise 



where the normalization constant Kq is determined by demanding that the total energy be 
equal to the injected energy: Kq = Eq/[Ex{2 + In {Ec/Ex))]- This spectrum is the same as 
that used by Protheroe, Stanev, & Berezinksy [10] and Jedamzik [12], and also agrees with 
the result of a detailed numerical integration of the full Bolzmann equation by Kawasaki and 
Moroi [9] It is a broken power law with a transition at E^ = Ex and a high-energy cutoff at 
E-y = Ec- We adopt the same energy limits as Kawasaki and Moroi, namely Ex = m^/SOT 
and Ec = ml/22T. Physically, these scales arise due to the competition between photon 
degradation rates. The scales rise as the temperature drops, in which case there exist more 
high-energy photons to break up nuclei. 

The zeroth-generation nonthermal photons then suffer additional interactions of Compton 
scattering, ordinary pair production off of nuclei, and 7 — 7 scattering. These slower pro¬ 
cesses further degrade the photon specturm. The evolution of the resulting “first-generation” 
photons is governed by 

= —Py(E,) - My(Ey)ry(Ey), (S) 

at Tx 

^For definiteness, we will refer to the decaying particle as a gravitino, though our analysis is general and 
pertains to any massive particle with electromagnetic decays in the lifetime range considered. 

^The spectrum stated in [9] includes further photon degradation, which has been factored out to determine 
our spectrum. 
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where nx = + 2 ;)^ exp (—t/rx) and tx are the decaying particle number density at 

redshift 2 : and mean lifetime, respectively. Also, is the photon energy spectrum, which is 
simply the product of the density of states and the occupation number fraction f^. Inte¬ 
grating ACy over all energies yields the number density of the injected photons. Further, 
r.y is the rate at which the photons are further degraded through further interactions with 
the background plasma. The key difference between and is that the rates degrading 
photons directly after injection are much faster than the rates that further degrade photon 
energy determining A/”^. We note that the effects due to the expansion of the universe on 
the photon spectrum are negligible because, during this epoch, electromagnetic interactions 
are much faster than the expansion rate. 

The dominant photon degradation rates are those for double photon scattering, Compton 
scattering and pair production off nuclei. Because their high rates are fast compared to 
the cosmic expansion, the photon distribution reaches quasi-static equilibrium (QSE). This 
distribution is given by setting (5) equal to zero, yielding 


- TO- 

This QSE solution is the same as that derived in [9], where it is called The photon 

spectrum can be determined easily from this equation, knowing that double-photon scat¬ 
tering dominates the high-energy region, whereas Compton scattering and pair production 
off nuclei dominate at lower energies. We recall that the redshift dependence of this QSE 
solution lies entirely in nx, Py, and T.^,. 


2.2 Photo-Destruction and -Production of Nuclei 


The equations governing the production and destruction of nuclei are very similar to those 
for photons, being given by 

^{Ea) = Ja{Ea) - Ua{Ea)Ta{Ea), (7) 

where Ja and T^ are the source and sink rates of primary species A. The derivative, d/dt 
takes into account the redshifting of energies and the dilution of particles due to the expansion 
of the universe. The source terms for the primary species are due to the photodissociation 
of background particles, and are dehned by: 


COO 

Ja{Ea) = / dE,M^^^{E,) a,+T^A{E,) 5 

rp J 0 




( 8 ) 


where Sa{E^) is the energy of the species produced by the photodissociation reaction 
7 -|- T A. The sinks are similarly defined by 


COO 

TAiEA) = E ldE^Af?^^iE^)(Ty+A^piE^). 


(9) 
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Since we are interested in calculating total abundances of elements, it is necessary to integrate 
(7) over the energy Ea- The equation then becomes 

^ = E a,+T^A{E,) -uaY. jdE,M^^^{E,) a,+A^p{E,). (10) 

This removes the redshifting term, leaving only the dilution term in the derivative, d/dt. It 
is useful to use the mole fraction 17 = ni/n-Q of baryons in a particular nuclide, rather than 
the absolute abundance. This allows us to take out the expansion effects, yielding 

= I dE,M^^^{E,) a,+T^A{E,) -YaY. jdE,M^^^{E,) a,pA^p{E,), ( 11 ) 

where np is the baryon number density and d/dt is an ordinary time derivative. It is also 
convenient to change from time differentiation to differentiation with respect to redshift, and 
also to extract the factor uxEq/tx out of In this way, we obtain 

dYA _ rMxv// / -1 \ 

dz HrTx V HrTx {I + z)"^ ) 

where we have used Eq = \Mx (corresponding to two-body decay into a photon and a particle 

of negligible mass), and have dehned r = n\/n/^ and Hr = ^327rGpr^^/3. We should note 
that, with these terms pulled out, the integrals become functions only of redshift, and the 
dependence on the properties of the decaying particle has been removed. This formulation 
is very useful in making the numerical implementation fast and efficient. 

2.3 Secondary Element Production 

In the previous Section we discussed the production of light elements by the photo-dissociation 
of heavier elements. However, the initial photo-production/destruction of light nuclei is not 
necessarily the only process that happens before thermalization. The primary interactions 
produce non-thermal particles, which then interact with the background plasma, degrading 
their energy. However, they may still have enough energy to initiate further, secondary 
nuclear interactions. 

We now modify the evolution of the primary particles described in the previous section 
to include energy-degrading interactions: 

^{Ea) = Ja{Ea) - Ma{Ea)Ta{Ea) - [hA{EA)MA{EA )], (15) 

at 


( 12 ) 

( 13 ) 

( 14 ) 
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Table 1: The relevant photo-dissociation reactions and their respective threshold energies are 
listed in the Table below, and their cross sections are listed in Appendix A. 


Reaction 

Threshold (E^^th) 

d( 7 ,n)p 

2.2246 MeV 

t( 7 ,n)d 

t( 7 ,np)n 

6.2572 MeV 
8.4818 MeV 

'^He( 7 ,p)d 

^He( 7 ,np)p 

5.4935 MeV 
7.7181 MeV 

^He(7,p)t 

^He(7,n)3He 

4He(7,d)d 

^He( 7 ,np)d 

19.8139 MeV 
20.5776 MeV 
23.8465 MeV 
26.0711 MeV 

®Li(7,np)'^He 

6Li(7,X)3A 

3.6989 MeV 
15.7947 MeV 

^Li(7,t)^He 

^Li(7,n)®Li 

’^Li(7,2np)^He 

2.4670 MeV 
7.2400 MeV 
10.9489 MeV 

^Be(7,^He)4He 

^Be( 7 ,p)®Li 

’^Be(7,2pn)'^He 

1.5866 MeV 
5.6058 MeV 
9.3047 MeV 
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where we have added the last term to include the energy degradation of the species A, where 
6 a = —dE/dt is the rate of energy loss. This term appears as an energy gradient, conserving 
particle number in the absence of sources and sinks. 

In most situations, the energy degradation rate is much faster than any sink, so that 
the sinks can be ignored. For unstable particles, if the lifetime of the particle is comparable 
to the stopping time of that species, then the Fa term cannot be ignored. In general the 
interactions are fast enough to reach a quasi-static equilibrium, but the form of the solution 
is somewhat more complicated than the photon case: 
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MK) 


ME)- 


( 16 ) 


Substituting Ja into this equation, we get 

poo 

Af^^^iEA) = , / _ dE^Af^{E^)a^+T^AiE^) exp 

Oa[-^a) t 


rSAiEj) 

Iea 


bA{EA)_ 


.(17) 


With this QSE solution in hand, we can determine the rate of secondary production of an 
element S: 


dMs 

dt 


{Es) = E pEAMf^{EA)TA+T'^s{EA)5 [Ss{Ea) - Eg ]. 

sec x' 


( 18 ) 


Integrating over Eg yields: 


dng 


dt 


sec 


poo 

dEAAfT^{EA)a^^^,^siEA)\vA\. 

. 7 J 0 

T 


( 19 ) 


Again, using the mole baryon fraction to remove expansion effects from the differential 
equation, we obtain 


dYg 


dt 


sec 


T.Yt' 

t' 


S(S.4)|W.4|. 


( 20 ) 


After some algebraic manipulation, we obtain the following evolution equation for secondary 
production: 


dYg 


dz 


' -( 1 - 1 - 2 :) exp 


-1 


I sec 
coo 


X / dE^ . ^ 

JszHea) kEqUx 


HrTx 
f Tx 


a^+T^A{E^) exp 


HrTx{l + Zf ) 

rSA(E^) 


Y, YtY^, / dE 


a 


A- 


A+T'^S 


iEA)\vA\ 


bA{EA) 


dE 


'Ea 


ta{e; 

^bAiE'A 


( 21 ) 


where r] is the baryon-to-photon ratio: f] = n^/n^. 

Table 2 lists the secondary reactions considered. Deuterium production does not occur 
within the lifetime range of interested to us, since the neutron decays before it has a chance to 
react with a proton to form deuterium, as pointed out in [7]. Also, we have verihed that the 
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Table 2: Relevant seeondary reaetions are listed below, in a format similar to Table 1. 


Reaction 

Threshold (Ep^th) 

p(n, 7 )d 

0.0000 MeV 

^He(t,n)^Li 

^He(3He,p)®Li 

8.3870 MeV 
7.0477 MeV 

4He(t,7)^Li 

^He(3He,7)^Be 

0.0000 MeV 
0.0000 MeV 


secondary production of mass-? elements is not significant, being small compared with that 
produced during BBN. The only signihcant secondary production is that of ®Li, as hrst shown 
in [12] and later in [13]. We show below that ®Li actually provides the strongest constraint 
for the lifetime range we are interested in. Note that the relevant threshold energies (Table 
2) are those in which the nuclei are in the cosmic rest frame, and thus are computed in the 
hxed-target laboratory frame. Consequently, these are about a factor of two higher than the 
center-of-mass thresholds used in [12] and [13]. 

3 Observational Constraints 

Before we discuss the results of our numerical analysis, we hrst discuss the current status of 
the observational determinations of the light-element abundances. The abundances subse¬ 
quent to any photo-destruction/production must ultimately be related to these observations. 
Furthermore, our results are dependent on the assumed baryon-to-photon ratio, which may 
either be determined through the concordance of the BBN-produced abundances or through 
the analysis of the CMB spectrum of anisotropies. As noted above, there is relatively good 
agreement between the two. 

3.1 Observed Light Element Abundances 

Through painstaking observations of very different astronomical environments, primordial 
abundances can be inferred for D, ^He, and ^Li. In addition, ^He and ®Li have also been 
measured, and can provide important supplementary constraints. Here we summarize the 
data and our adopted limits: more detailed reviews appear in [15]. For all nuclides, accurate 
abundance measurements are challenging to obtain, due to systematic effects which arise 
from, e.g., an imperfect understanding of the astrophysical settings in which the observations 
are made, and from the process by which an abundance is inferred from an observed line 
strength. 

Deuterium is measured in high-redshift QSO absorption line systems via its isotopic shift 
from hydrogen. In several absorbers of moderate column density (Lyman-limit systems), D 




has been observed in multiple Lyman transitions [42, 43]. Restricting our attention to the 
three most reliable regions [42], we hnd a weighted mean of 

= (2.9±0.3) X 10"l (22) 

VH/p 

We note, however, that the per degree of freedom is rather poor ( ~ 3.4), and that the 
unweighted dispersion of these data is ~ 0.6 x 10“^. This already points to the dominance 
of systematic effects. Observation of D in systems with higher column density (damped 
systems) hnd lower D/H [44], at a level inconsistent with (22), further suggesting that 
systematic effects dominate the error budget [45]. If we used all hve available observations, 
we would hnd D/H = (2.6 ±0.3) x 10“^ with an even worse per degree of freedom (~ 4.3) 
and an unweighted dispersion of 0.8. As an upper limit to D/H, we adopt the 2-(j upper 
limit to the highest D/H value reliably observed, which is D/H = (4.0 ± 0.65) x 10“®, since 
we cannot dehnitively exclude the possibility that some D/H destruction has occurred in the 
other systems. 

We also require a lower limit on the primordial D abundance. Since Galactic processes 
only destroy D, its present abundance in the interstellar medium[46], D/H = (1.5±0.1) x 10“^ 
provides an extreme lower limit on the primordial value, which is consistent with (22). 
Therefore, we adopt the limits 

1.3x10-® <5.3x10-®. (23) 

This lower bound is quite conservative, in light of the fact that the existence of heavy elements 
conhrms that stellar processing and thus D destruction has certainly occurred at some level. 

Unlike D, ^He is made in stars, and thus co-produced with heavy elements. Hence the 
best sites for determining the primordial ^He abundance are in metal-poor regions of hot, 
ionized gas in nearby external galaxies (extragalactic HII regions). Helium indeed shows a 
linear correlation with metallicity in these systems, and the extrapolation to zero metallicity 
gives the primordial abundance (baryonic mass fraction) [47] 

Yp = 0.238 ± 0.002 ± 0.005. (24) 

Here, the hrst error is statistical and reflects the large sample of systems, whilst the second 
error is systematic and dominates. 

The systematic uncertainties in these observations have not been thoroughly explored to 
date [48]. In particular, there may be reason to suspect that the above primordial abundance 
will be increased due to effects such as underlying stellar absorption in the HII regions. We 
note that other analyses give similar results: Yp = 0.244 ± 0.002 ± 0.005 [49] and 0.239 
±0.002 [50]. For concreteness, we use the ^He abundance in (24) to obtain the range 

0.227 <Yp< 0.249, (25) 

taking the 2-a range with errors added in quadrature. 

Helium-3 can be measured through its hyperhne emission in the radio band, and has 
been observed in HII regions in our Galaxy. These observations hnd [51] that there are no 
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obvious trends in ^He with metallicity and location in the Galaxy, but rather a ^He ‘plateau’. 
There is, however, considerable scatter in the data by a factor ~ 2, some of which may be 
real. Unfortunately, the stellar and Galactic evolution of ^He is not yet sufficiently well 
understood to conhrm whether ^He is increasing or decreasing from its primordial value [23]. 
Gonsequently, it is unclear whether the observed ^He ‘plateau’ (if it is such) represents an 
upper or lower limit to the primordial value. Therefore, we do not use ^He abundance as 
a constraint. If future observations of ^He could hrmly establish the nature of its Galactic 
evolution, then ^He could be restored as a useful constraint on decaying particles, particularly 
in concert with D. 

The primordial ^Li abundance comes from measurements in the atmospheres of primitive 
(Population II) stars in the stellar halo of our Galaxy. The ^Li/H abundance is found to 
be constant for stars with low metallicity, indicating a primordial component, and a recent 
determination gives 


= (1.23 ± 0.06l[];|) X 10“^° (95% GL), (26) 

p 

where the small statistical error is overshadowed by systematic uncertainties [52]. The range 
(26) may, however, be underestimated, as a recent determination [53] uses a different proce¬ 
dure to determine stellar atmosphere parameters, and gives ^Li/Hp = (2.19 ± 0.28) x 10“^*’. 
At this stage, it is not possible to determine which method of analysis is more accurate, 
indicating the likelihood that the upper systematic uncertainty in (26) has been underesti¬ 
mated. Thus, in order to obtain a conservative bound from ^Li, we take the lower bound 
(once again combining the statistical and systematic errors in quadrature) from (26) and the 
upper bound from [53], giving 



9.0x10“^^ < 



< 2.8x10“^°. 

p 


(27) 


Finally, ®Li is also measured in halo stars, in which the ®Li/^Li ratio is inferred from 
the (thermally blended) isotopic line splitting. The lowest ®Li abudances comes from stars 
with primordial Li, which yield ®Li/^Li = 0.05 ± 0.01 [54]. The ®Li in these stars is not 
primordial, as it is produced by cosmic-ray interactions with the interstellar medium [55], 
predominantly aa ®Li -!-■■■. These same processes lead to the production of Be and B, 
which are observed in halo stars at levels consistent with ®Li cosmic-ray production. Since 
the observed ®Li abundances are consistent with being entirely Galactic in origin, we can 
use these to set an extreme upper limit on the primordial ®Li abundance. One complication 
enters, due to the smaller binding energy of ®Li relative to ^Li. This means that ®Li could in 
principle suffer depletion in stars due to nuclear burning, without a similar depletion of ^Li. 
However, once nuclear burning becomes effective, ®Li depletion factors become extrelemly 
large making such observations extremely unlikely [56]. It is therefore safe to use the 2-a 
upper bound on the ®Li/^Li ratio 


( 



< 0.07. 


(28) 
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Other depletion processes such as diffusion (included in the estimate of systematic uncer¬ 
tainties in (26)), would affect both ®Li and ^Li similarly and not their ratio. It is also useful 
to consider the upper bound on ®Li/H alone 

<2xl0-ih (29) 

3.2 Cosmic Microwave Background Anisotropy Measurements 

Cosmic Microwave Background (CMB) anisotropy data are now reaching the precision where 
they can provide an accurate measure of the cosmic baryon content. Given a CMB measure¬ 
ment of rj, one can use BBN to make dehnite predictions of the light element abundances, 
which can then be compared with the observations discussed above. This comparison con¬ 
strains the effects of decaying particles more powerfully than if only the BBN calculations 
were available to constrain rj. 

Recent results from DASI [18] and CBI [19] indicate that = 0.022to;oo3, while 

BOOMERanG-98 [20] gives = 0.02llo;oo3- These determinations are somewhat lower 
than the central values found by MAXIMA-1 [21]: = 0.026to;oo6 and VSA [22]: = 


0.029 ± 0.009. Taking a CMB value of 

= 0.022 ± 0.003 or ? 7 io,cmb = 6.0 ± 0.8 (30) 

at the 1-a level, we would predict the following light element abundances: 

^He : 0.248 ± 0.001 (68% CL) (31) 

D/H X 10^ : 2.7toi (68% CL) (32) 

^He/H X 10^ : 0.9 ± 0.1 (68% CL) (33) 

^Li/H X 10^° : 3.4;j| (68% CL) (34) 

Note that these numbers are not outputs of BBN calculations corresponding to rjiQ = 6.0, 


but rather are the peak values of a likelihood function found by convolving the results of 
the BBN Monte Carlo with an assumed Gaussian for the distribution of CMB rj values. For 
further details, see [16, 17]. With MAP data, the accuracy of rjcmh should be 10% or better, 
which will give even tighter predictions on the light elements. 

4 Model Results 

We have implemented numerically the decaying-particle cascades discussed in Section 2. 
Using BBN light-element abundance predictions [16] as initial conditions, we calculate the 
hnal abundances for particular sets of baryon and dark matter parameters. The three free 
parameters are: 

Cx = — wdTx = r Mx = 2r Eq, (35) 

Tx and rj. 
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4.1 Analytic Discussion 

Some simple analytic approximations allow us to gain insight into the essential physics in 
our problem. As we will see, the following analytic treatment reproduces well the behavior 
of the light element abundance mountains and deserts in our parameter space. 

The dependence on tx can be understood [11] in terms of the characteristic energy scales 
in the photon spectrum (4). Both the break Ex and the cutoff Ec scale as Ei oc 1/T. 
Thus, in the ‘uniform decay’ approximation where all particles decay at t = tx, the decay 
occurs at T ~ 10“"^MeV (rx/10® s)“^/^. Consequently, we have Ex ~ 28MeV (rx/10® s)^/^, 
and Ec ~ 103 MeV (rx/10® s)^/^, and cutoffs thus increase with tx- In other words, higher- 
energy photo-erosion processes can occur for longer lifetime values. Comparison with Table 1 
shows that, as tx increases, first Eq and then Ex pass the threshold energy for any 
given process, at which point the process becomes important. A reaction can turn on when 
Ec ^ Eth, which in the uniform decay approximation occurs when 


Tx > 10® S 


E, 


th 


10 MeV 


(36) 


while for shorter tx the channel is closed. The reaction grows stronger when Ex ^ 
which occurs when 


Tx ~ lO’^ S 


Eth 

10 MeV 


2 


(37) 


We can also understand the (x dependence of photodestruction (and secondary production) 
analytically, as follows. In the limit of small (x, the decaying particle has no influence on the 
light-element abundances as predicted by primordial nucleosynthesis, predicting a universe 
made of mostly hydrogen and ^He, with small but signihcant amounts of D, ^He, and ^Li. 
Lithium-6 is not produced in signihcant quantities. Going beyond this trivial case we use a 
similar treatment as above, and employ the uniform decay approximation. To begin, as long 
as a reaction can proceed, a typical shower photon has energy 


(E) = 56 MeV 


Eth \ 

10 MeV; 



(38) 


SO that the number of such photons per decay is ~ T^o/ {E). Had the lower-energy piece of 
the power law been much steeper (i.e., with a power index p > 2) we would have (E) ~ Eth, 
and if it was much shallower (i.e., p < 1) we would have (E) ~ Ex- However, we he in an 
interesting regime where (E) ~ {Eth/Ex)^~^Ex, where 1 < p < 2. Thus the nonthermal 
photon density is 


= N^nx 


Cxn^^ 

2{E) ' 


(39) 


These photons are thermalized at a rate per photon of Ttherm ~ ne<JT ~ nB<TT- The rate 
per photon for the photodestruction of species T to yield species P is T 

T^P ^ Tlx^T^P — 

therm, where Yt = TiT/n-B, and fr^p = ap^p/cp is the relative strength of the cross 


12 



section for photodestruction of T into P, compared with the thermalization cross section, 
which we take as the Thompson cross section for this discussion. Consequently, the change 
in the number density of A is the net production rate per volume (Tr^^ — rA^p)n'^^ times 
the loss time Th^erm) 


duA = uboYa ~ ---. 

-L therm 

We see that the fractional change in an abundance is given by 

Cx(Yt. . 

Ya ^ 2r]{E) [ya^ 



(40) 


(41) 


If we look at the two extremes when either production or destruction dominates, we can 
derive the behavior of (x, given that the fractional change in the abundance is 1/2: 


^prod 


^dest 


rsj 


rsj 


3.2x10-^^ GeV 



6.3x10“® 


/S.OxlO-^ 


5.0x10 /?7io\ / Eth y 

fp^A )\q)\ 20 Mev) 

^ Y ( 

/ V 6 / V2.224 MeV; UO^ s 


TX Y 

108 S/ ’ 

1 

4 


( 42 ) 

( 43 ) 


where the numbers are those appropriate for D. 

This same treatment can be extended to secondary production of light elements. Since ®Li 
is the only signihcant secondary production, this is the only example we consider here. In this 
case, only a fraction of the primary products have enough energy for this reaction to proceed, 
because of interactions with the background plasma. We estimate as follows the fraction that 
can react to form ®Li. Each prospective reactant is produced with initial energy Sa{{E)), and 
the total amount of energy that can be lost between collisions is Sioss = bA{{E))/nBaT\vA\- 
The fraction of reactants left is the ratio of these two energies, Np ~ Sa{{E))/S ioss ~ 0.001. 
Thus the number density of these remaining non-thermal particles is 


n 


inj 

p 


NpSuA 


2ri{E) 


(44) 


The rate per particle can be described in a similar way as before, where T^^s = fA^sYp'TthermP- 
Since these particles are non-relativistic, there is a factor P = v/c in the interaction rate. 
We can thus determine the change in abundance of the secondary species: 


Sns ~ nsSYs ~ Up • 

therm 

The fractional change in the abundance is then given by 

^>5 




^ CxYtYt' f f 
2ri{E)Ys 


(45) 


(46) 
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Figure 1: Contours of the ^He mass fraction Yp (a) in the plane, for rjiQ = 6, and 

(b) in the {Cx,v) plane, for tx = 10® sec. See the text for discussion. 


Using parameters appropriate for ®Li and fd ~ 0.01, we derive the value of C,x when secondary 
production becomes important: 


~ 2.6x10“^^ GeV 


2.5x10^^^ 


5.0x10 / 0.05 \ /T]io\ / Eth y \ 

fr^A ) [fA^sJ V 6 A 20 MeV; Uo® s) 


.(47) 


One should note that, at the high photon energies required to induce these tertiary reac¬ 
tions, double-photon scattering is comparable to the Compton scattering and nuclear pair- 
production mechanisms for photon energy-loss. This weakens the dependence of the reaction 
rates on the baryon density. 

We now turn to the full numerical results, using these analytical estimates as a guide to 
interpretation. 


4.2 Numerical Results 

Our numerical results come from the integration of (12) and (21), using the input spectra 
of (4). We present our results by showing abundance contours in (x — Tv space at hxed 
r/io = 6, and in (x — V space at hxed tx = 10® s. A key innovation of the present work is a 
detailed htting of the energy dependence of the relevant cross sections. These are discussed 
in Appendix A, which also contains htting formulae. 
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Figure 2: Deuterium abundance contours plotted as in Fig. 1. 


We begin our discussion with the most abundant compound nucleus, namely “^He. Since 
the other light elements are predicted by BBN to have abundances that are orders of magni¬ 
tude smaller than ^He, no signihcant production of ^He can take place. Thus electromagnetic 
showers from decaying particles can only destroy ^He. These photodestruction processes have 
an energy threshold of F^th ~ 20 MeV and so, from (36), we expect this process to become in¬ 
efficient for Tx ^ 4 X 10® s, and to shut down completely when rx ^ 4 x 10® seconds. Indeed, 
this is what is seen in Fig. 1(a), where we plot contours of Yp in the (Cx, Fv) plane for r/io = 6, 
the value preferred by CMB analyses (30). For rx ^ 4 x 10® s, the ^He destruction factor 
goes from a small perturbation to a large one as C,x grows from 10“^® to 10“® GeV, until the 
region Qx ^ 10“® GeV becomes a "^He ‘desert’. Over the region rx ~ 4 x 10® s — 4 x 10® s, 
'^He destruction becomes important only at increasingly high (x- This general behavior has 
an impact on all of the other light elements, as ^He is the only important source for them. 

In Fig. 1(b), we again plot contours of Yp, but now in the ((x,p) plane for tx = 10® s. 
We see the generic features mentioned above, that for low (x ^He is at its BBN predicted 
value and for very large (x it is destroyed. At large ?), the non-thermal photons are more 
quickly thermalized, thus having less energy on average and making it more difficult to 
destroy nuclei. This linear rise is as predicted in (42). 

Fig. 2 plots the corresponding contours of D/H. Several distinct regions are apparent. For 
all Tx, in the region (x ^ 10“^^ GeV, the decaying particles make only a small perturbation 
to the primordial value of D/H. For values of (x ^ 10“^^ GeV, the decaying particles can 
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lead to significant perturbations to D/H, the sign and magnitude of which depend strongly 
on Tx- In the case of D, the important production processes are due to "^He + 7 : the other 
light elements have negligible abundances compared to "^He, and thus are unimportant as 
D sources. The D production channels have thresholds i?th ~ 25 MeV, so we expect these 
process to become inefficient for tx ^ 5 x 10® s. Furthermore, production is only efficient 
when ^He destruction is not large, i.e., when ^ 10“® GeV. Of course, D production can 
only occur when there is sufficient ^He destruction, so we expect signihcant production to 
occur only up to some maximum Qx- For higher C,x, photo-destruction becomes so dominant 
that any production yields are immediately broken up via further photo-destruction reac¬ 
tions, leaving a universe hlled with only protons. Thus, as C,x increases, the D/H abundance 
rapidly declines, dropping to and then below its primordial abundance, to approach zero in 
a D ‘desert’. Thus, we expect D production only in a region bounded from above and below 
in Cx, and to the left by rx ^ 4 x 10® s. These expectations are met by the D/H ‘mountain’ 
in Fig. 2, which stretches between Cx = 10“^® GeV to Cx = 10“® GeV. 

For Tx ^ 5 X 10® s, D production from ^He does not occur. Since the photo-destruction 
processes 1 ^( 7 , n)p has a threshold of i?th = 2.224 MeV (the D binding energy), D destruction 
drops out at Tx ^ 4 X 10® s, as seen in Fig. 2. Finally, note that the competing processes of 
D production and destruction balance for some regions of parameter space, where D retains 
its primordial values. These are the ‘channels’ which separate the regions we have already 
discussed. 

In Fig. 2(b), the contours of D/H are shown in the (Cx, 7 ) plane for tx = 10® s. Just like 
"‘He, these D contours show the features sketched out in our analytic discussion. For low Cx, D 
is at its BBN value. At intermediate values, we are climbing the photo-production mountain. 
At even higher values of C,x, we enter the photo-erosion desert. Again the dependences on rj 
are entirely due to the photon energy-loss mechanism being more efficient at higher baryon 
density. 

Results for ^Li/H are shown in Fig. 3. We see that the results are qualitatively similar to 
those for ^He, reflecting the fact that secondary ^Li production is negligible. Since ^Li is only 
destroyed, its weak binding compared to ^He leads to the wider expanse of the ^Li ‘desert’. 
As mentioned earlier, we considered the secondary production of ^Li, through the reactions 
^He(t, 7 )^Li and ^He(®He, 7 )^Be. These reactions have no energy threshold, only a strong 
Goulomb barrier, so a priori they would seem important. However, the net production of 
^Li through these reactions is small compared to its primordial value set by BBN. This is 
due to the small cross sections for mass-7 production. 

Unlike ^Li, standard BBN does not produce ®Li in any observable quantity, so any other 
production mechanism is important. The behavior of ®Li, seen in Fig. 4a, can be understood 
in terms of the ^He and ^Li dynamics, because these are the two ®Li sources. As with 
D, we see a roughly horizontal ‘mountain’ of ®Li production, which is bounded on the left 
by threshold effects. The dominant production channel is from secondary reactions, and 
thus is tied to the "^He destruction threshold. The secondary reactions ^He(t,n)®Li and 
^He(®He,p)®Li, also have low energy resonances, further increasing their yields. The ®Li 
production from ^Li and ^Be, however, has lower thresholds, and thus becomes dominant for 
small lifetimes, tx ^ 10^ s. These channels are important where ^Li destruction is moderate 
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Figure 3: Contours of'^Li/H, plotted as in Fig. 1. 
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Figure 4: Contours of^hi/H, plotted as in Fig. 1. 
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Figure 5: Contours of ^Li/^Li, plotted as in Fig. 1. 


but still sufficient to make significant amounts of ®Li, leading to a break in the slope of the 
®Li curve at C,x ~ 10“^^. The inclusion of two ®Li destruction channels leads to the ‘nose’ 
at Cx r\j 10 ^ GeV and tx ~ 10^ s. This feature is absent in the ®Li plot in [13], since there 
only a single destruction channel was considered. 

In Fig. 4b we see that the secondary production dominates the evolution of ®Li for quite 
low C,x- Since the ®Li secondary-production cross sections have thresholds, the average 
initiating photon energy must be higher than in the standard photo-destruction process. 
In this higher-energy regime, double-photon scattering is comparable to the other photon 
energy-loss mechanisms. This reduces the dependence on the baryon density, effectively 
flattening out the contours in the {Cx,v) plane where secondary production is important. 
At higher Cx, photo-destruction of ®Li takes over, taking us to the ®Li desert. 

As noted in [11], the ®Li/^Li ratio offers additional constraints besides those provided 
by each nuclide separately. We plot ®Li/^Li contours in Fig. 5, where we see that these are 
smoother than the contours for ®Li alone. Because ®Li dissociation has a higher threshold 
than ^Li - see Table 1 - ®Li/^Li is large in the ‘desert’ region of high (x, though each 
individual abundance is quite small. For smaller (x, ^Li increases with (x, while ^Li either 
remains constant or decreases. The upshot is that ®Li/’^Li grows with (x, as seen in Fig. 5. 

Having described the physics that leads to the abundance patterns we have computed, we 
can now discuss how this physics allows the observed abundances to place constraints on de¬ 
caying particles. To obtain these constraints, we combine this analysis with the observational 
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data discussed previously. 


5 Limits on Unstable Relic Particles 

We now impose the observed light-element abundance constraints of Section 3.1 with the 
results of the previous section. We do this for each element individually, then combine the 
results to obtain the strongest constraints. 

We remind the reader that light element constraints on decaying particles depend on 
7], but of course one cannot use the standard BBN limits on rj as part of ones limits. In 
the past, this difficulty has only been overcome by adopting limits on rj derived from non- 
nucleosynthetic arguments. These limits have, until recently, been rather weak, which has 
weakened the power of the light element constraints. This situation has now changed dras¬ 
tically. We recall that the CMB now imposes rjiQ ~ 6 with rather small uncertainty. Thus, 
if we adopt the CMB results, we no longer must treat rj as weakly constrained by non-BBN 
arguments, strengthening the results we derive. 

We quote limits for rjiQ = 6, and emphasize results for tx = 10® s, which is roughly the 
lifetime for which the constraints are the strongest, and is also within the range of current 
interest for gravitino decays. Results at lower tx weaken rapidly, while the constraints at 
higher tx scale roughly as as in (42), (43) and (47). 

Also, the current observational status of standard BBN comes into play. Namely, the 
present observational data on ^He and ^Li are in tension with those for D, the former pre¬ 
ferring ? 7 io ~ 3 and the latter rjiQ ~ 6. On the one hand, if standard BBN is correct and 
the tension is due to systematic errors, the result is that these errors weaken the constraints 
that one can place on decaying particles. On the other hand, if the observations were to 
improve to the point that the light-element disagreement can no longer be accommodated, 
this could herald new physics. In this case, decaying particles offer a way [11] of reconciling 
the abundances of the light elements, in which case one also derives estimates of the required 
Cx and Tx- 

We turn hrst to the elements which are only destroyed, namely ^He and ’^Li. The observed 
constraints on ^He (25) give 

Cx(^He) < 2.5 X 10"^° GeV, (48) 

which is driven by the lower limit Yp > 0.227 that we adopted. The sharp drop of ^He with 
increasing Qx (he., the descent into the desert of Fig. 1) ensures that the constraint on (^x 
is insensitive to the precise Yp limit chosen. The situation is similar for ^Li, for which (27) 
give the weaker constraint 

Cx(^Li) < 2x10“® GeV. (49) 

For deuterium, net production and net destruction are both possible. In terms of Fig. 
2, this means that limits on D exclude the ridge in the D mountain, while allowing regions 
at higher and lower C,x- In particular, the observed D abundances (22) allow the range 
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Figure 6: Exclusion limits based on deuterium are shown in (a) the (Cx, tx) plane forrjiQ = 6 
(b) the {Cx,v) plane for tx = 10® s. The dark (blue) shaded region corresponds to an 
overabundance of D/H, while the light (blue) shaded region corresponds to an underabundance 
ofD/H. 

2x10“® GeV < Cx ^ 3x10“® GeV, but the "^He and ^Li constraints are each able to exclude 
this regime. Gonsequently, the only remaining region is the low-(' side of the mountain, 

Cx(D) < 3x10“^^ GeV. (50) 

The since ®Li is not produced signihcantly in standard BBN, only production is important 
for low (x, while for higher (x destruction dominates. Thus, the situation is similar to that 
of D: there is a ®Li mountain, with the observations allowing a narrow high-(Cx region and a 
large low-^x region. The ®Li/H abundance of (29) gives 

Cx(®Li) < 5x10“^^ GeV (51) 

in addition to a higher region that is discordant with ^He and ^Li. The ®Li/^Li ratio (28) 
gives 

Cx(®Li/^Li) < 7.0x10“^^ GeV. (52) 

Figure 6 summarizes our results for the constraints based on D/H in both the (Cxj'T’x) 
plane (for r/io = 6) and the {Cx,v) plane for tx = 10® s. The dark (blue) shaded regions 
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Figure 7: As in Figure 6, including the constraints from “^He - medium (pink) shading and 
^Li - medium-light (green) shading. 


correspond to an overabundance of D/H, i.e. regions where there is net production of deu¬ 
terium. The lighter (blue) shaded regions represent an underabundance of D/H or regions 
where there is net destruction of deuterium. Notice that the thin strips which are unshaded 
for which the D/H abundance is acceptable. These will be excluded when the constraints 
from the other light elements are included. 

In Figure 7, we include the constraints from ^He and ^Li. Here we superimpose the “^He 
constraint, shown as the medium shaded (pink) region, and the ^Li constraint, shown as the 
medium-light (green) shaded region, on the D/H constraints. We see that D, ^He, and ^Li 
alone, i.e., primordial species, impose a limit of 

Cx <Cmax = 3.5 X 10-^1 GeV, (53) 

which is dominated by the limits from D. We can do better if we include ®Li. Our limit 
®Li/H < 2 X 10“^^ pushes the above constraint down to 

CX <Cmax = 5x10-12 Qgy ^ 54 ) 

for Tx = 10^ s as seen in Figure 8 by the dark (red) shaded region. The constraint from 
the ®Li/i^Li ratio is shown as the light (yellow) shaded region. Notice that it becomes the 
stronger constraint at rj < 5.0. 
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Figure 8: ^45 in Figures 6 and 1 including the constraints from ®Li (dark red). 


These constraints are subject to uncertainties in the ®Li limit, due both to the possible 
stellar depletion of ®Li and the known Galactic production of ®Li by cosmic rays. Our limit is 
intended conservatively to allow for both effects. Even so, we see the power of ®Li. We thus 
urge further observations of the Li isotopic ratio, as a hrmer understanding of this nuclide 
could further strengthen the constraint we have derived. 

As already noted, if the observed light-element abundances retain their current central 
values, but the error budget shrinks, then the light-element data will be in discord with 
standard BBN. Decaying particles might provide one possible means of reconciling such 
light-element observations and theory. As an illustrative example, consider the case in which 
the CMB fixes r)io = 6, and the observed light element abundances remain as above, but 
with the total error budget equal to that of the current statistical errors. Then ^Li and D 
would be in signihcant disagreement. One could, however, bring these nuclides and ^He into 
agreement by appealing to the decaying-particle scenario we have laid out here, the allowed 
region of parameter space still open being one in which a non-zero C,x is preferred. The 
new ^Li upper limit would eliminate low values of allowing only a narrow band with 
C 10“®. The observations would force us to live in the narrow channel where D production 
and destruction are nearly balanced, with a decaying particle lifetime rx ~ 3 x 10® sec. In 
this regime, ^He is at its BBN value, because of its high photoerosion threshold, shown 
in Table 1. However, the more weakly bound ^Li is destroyed, at just the right level to 
bring the observations in accord with the D observations. This rather fine-tuned scenario is 
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testable with ®Li observations. In this region of parameter space, we predict a ®Li abundance 
of ®Li/H ^ 10“^^ and ®Li/^Li ^ 0.03. This ®Li abundance would appear as a pre-Galactic 
plateau in halo stars. Indeed, the ®Li level would be large enough to dominate the cosmic-ray 
component of ®Li over most of the Population II metallicity range. 


6 Application to Cosmological Gravitinos 


We now illustrate the impact of our calculations by discussing their implications for cosmo¬ 
logical gravitinos. In conventional supergravity scenarios, the gravitino is expected to have 
a mass comparable to that of supersymmetric partners of standard Model particles, which 
should weigh less than about 1 TeV if they are to stabilize the gauge hierarchy [57]. There¬ 
fore, the gravitino is usually thought to weigh between about 100 GeV to 10 TeV, though 
both larger and smaller masses have sometimes been considered. The lightest supersymmet¬ 
ric particle (LSP) is generally thought to be the lightest neutralino y, a model-dependent 
mixture of the photino 7 , the zino Z and the neutral Higgsinos H\ 2 [58]. The LSP would 
be stable in models in which R parity is conserved, as we assume here. On the other hand, 
the gravitino would be unstable, with a partial G —y -|- 7 decay rate calculated to be 

1 (M%\ , 

r(G - X7) = 4 Ol,, (55) 


where Oy^.y is the fraction of 7 in the wave function of the LSP y, and Mp is the Planck 
mass. In many models, the LSP is essentially a pure U{1) gaugino (Bino) B, in which case 
= cos"^ 9w and the neutralino mixing does not suppress the gravitino decay rate (55). 
We assume for now that y = i? and that no other gravitino decay modes are signihcant, in 
which case the gravitino lifetime is 




2.9 X 10® 


100 GeV\^ 


(56) 


We discuss later the modihcations to our analysis needed if the LSP is essentially a pure 
Higgsino, another possibility sometimes considered, or if other decay modes are open to the 
gravitino. 

The production of gravitinos in the early Universe has been the subject of heated discus¬ 
sion. An unavoidable contribution is thermal production [59, 8 , 60, 61]. Here, we will apply 
our results in combination with the recent calculation in [61], which gives 


y- = 


1.2 X 10"^^ 



X 


Ti 


R 


1010 GeV 


(57) 


where Tr is the maximum temperature reached in the early Universe^. In conventional in¬ 
flationary cosmology, Tr is the reheating temperature achieved at the end of the inflationary 

^We note that this calculation is based on the dominant strong contributions to gravitino production. 
Electroweak corrections would enhance the production rate by about 5 - 20 %. 
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epoch. In some inflationary scenarios, there may be additional gravitino production, either 
during the inflationary epoch or later, before thermalization is achieved. Either of these ef¬ 
fects would only accentuate the potential problem we discuss below, and we do not consider 
such possibilities here. 

In (57), m~ is the low-energy gluino mass. In supersymmetric models with gaugino mass 
unification, there is a definite relation between the gluino mass and the bino mass, which, at 
the one-loop order sufficient for our purposes, is m~/m^ ~ Q; 3 (m~)/Q;i. Typically, this ratio 
is between 5 and 6. If the bino is the LSP, as we are presently considering, then ttIq > mg. 
The middle term in (57) is therefore never larger than 4, and tends to unity for large gravitino 
masses. Thus we estimate 


Eg ~ (0.7-2.7) X 10"^^ X 


I 

101» GeV/ 


(58) 


in this case. 

The upper limit (1) can be expressed as a limit on Yg: 


Eg < 5 X 10"^^ X 


/ lOO GeV 

V 


(59) 


Comparing the calculated abundance (58) with the upper limit (59), we infer the following 
upper limit on the reheating temperature Tr, for Mg ~ 100 GeV: 

Tr < (1.9 - 7.5) X 10^ GeV (60) 

This upper limit is far smaller than the reheating temperature Tr ~ 10^^ GeV expected in 
conventional inflationary scenarios. As noted in the introduction, this bound places impor¬ 
tant constraints on models of baryo/leptogenesis. 

We recall that the upper limit (60) comes from a combination of data and calculations of 
different light-element abundances. It could not be obtained by considering the deuterium 
abundance alone, as this would allow a ‘tail’ of the parameter space extending to large 
gravitino abundances Yg, as well as an isolated ‘channel’ at large Yg. The ‘tail’ cannot 
be excluded just by considering also the abundance of ^He, though the ‘channel’ probably 
can. However, as discussed earlier, the ‘tail’ can be excluded by also considering the ®Li 
abundance. We have already discussed why we think that measurements of the relevant 
photoreaction cross sections and the astrophysical data are now sufficiently reliable for the 
®Li data to be regarded as a serious constraint. 

As we commented at the end of the previous Section, if the observed light-element abun¬ 
dances were to retain their current central values while the error budget shrank, the light- 
element data would become inconsistent with standard BBN, and decaying particles might 
be able to reconcile such light-element observations with theory. This scenario would predict 
a lifetime rg ~ 3x10® sec and an abundance Cg ~ 10“® GeV. These constraints would place 
the gravitino mass at Mg 460 GeV, with a relic abundance of Yg ^ 2x10“^^, corresponding 
in turn to a reheating temperature Tr ~ (0.8 — 3.1) x 10® GeV. 

How might the potentially embarrassing eonclusion (60) he avoided or evaded? 
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The first option one might consider is diluting the density of gravitinos by several orders 
of magnitude some time between their production at a temperature close to Tr and the 
period when they decay. This large entropy release should certainly occur before BBN, i.e., 
when the age f ^ 1 s, in order to avoid destroying its predictions completely. Very likely, such 
a large entropy release would also have had to occur before baryogenesis, in order to avoid an 
unacceptable dilution of the primordially-generated baryon asymmetry. The latest epoch at 
which baryogenesis is seems likely to have occurred is the electroweak phase transition, which 
occurred when the age t ~ 10“^° s. Afleck-Dine baryogenesis [40] offers one such possibility. 
In these models the Universe becomes dominated by the oscillation of a scalar field along 
a supersymmetric flat direction. In general, the net baryon asymmetry produced in these 
models can be quite large, actually necessitating late entropy production [62]. The dilution 
of the gravitino abundance would be an immediate consequence. The issue of gravitino 
production and dilution in connection with baryogenesis was recently considered [63] in the 
context of the pre-big bang scenario [64]. 

It is possible that the mass and decay modes of the gravitino are such that the bounds 
discussed here, which cover lifetimes from 10"^ - 10^^ s, become inapplicable. We see from 
(55) that the rate for G —X -f 7 decay would increase by three orders of magnitude if 
Mq were one order of magnitude larger. In fact, a heavier gravitino might have additional 
decay modes open kinematically, possibly decreasing its lifetime by another two orders of 
magnitude if all the MSSM particle weighed less than Mq. If Tq ~ 10^ s, as might occur 
if Mq ~ 1 TeV and it could decay into Xew + Xew as well as x + 7 , where Xew is 
any Standard Model particle with only electroweak interactions, our limits are signihcantly 
weakened. However, we remind the reader that outside the range 10^ - 10^^ s other bounds 
come into play. At smaller lifetimes, hadronic decay products will upset the prediction of 
BBN [5, 6 , 14]. From the recent results of [14], one hnds that for lifetimes in the range 1 - 
10^ s, the upper limit to Yq is of order 10“^^, and this bound weakens at lower lifetimes. For 
r < few X 10“^ s, the BBN limit disappears. Similarly, at longer lifetimes > 10^^ s, there 
are non-negligible constraints from the observed gamma-ray background [65, 7, 66 ]. These 
are strongest for decay lifetimes of order the age of the Universe ~ 10^® s, where the bound 
on ttIqYq is of order 10 “^®. At still longer lifetimes, the bound weaks quickly and becomes 
ineffective at lifetimes longer than about 10 ^^ s. 

Could the cosmological embarrassment be avoided if Mq -C 100 GeV? In such a case, the 
gravitino would presumably be the LSP, and absolutely stable if R parity is conserved, as we 
have been assuming. In fact, as argued in [61], a reheating temperature of order 10^'^ GeV 
would result in an acceptably large relic density of gravitinos. Converting the gravitino 
abundance in (58) to its contribution to closure density, one hnds. 




(0.026 


0.1) X 


( \ 

VlOOGeV; 


10‘» GeV) 


(61) 


In this case, the lightest MSSM sparticle (presumably the lightest neutralino x) would be 
the next-to-lightest supersymmetric particle (NLSP), and would itself be unstable: the decay 
NLSP — 7 -|- G would have a lifetime similar to (56), and the bound (57) could be applied 


26 



( 62 ) 


to the NLSP abundance Y^lsp, which in terms of is 
< 2 X 10"^ 

However, suppressing to such a low value seems very difficult. Characteristic values of 
in the MSSM are (9(0.01—10.0) [67], corresponding to ~ 10“^°, which is far above the 
bound (57). Indeed, in much of the phenomenologically allowed supersymmetric parameter 
space, the relic density is too large. Relic neutralino densities as low as (62) would require 
special parameter choices for which the neutralino is either primarily a massive Higgsino, in 
which case important coannihilation effects might suppress the abundance, or the neutralino 
mass happens to be very close to half the mass of the Heavy Higgs scalar and pseudo scalar, 
so that rapid s-channel annihilation reduces the neutralino density. In the former case, it is 
in fact unlikely that coannihilations are sufficiently strong to satisfy the bound (62), despite 
reducing the relic density to a level well below critical. In the latter case, however, very low 
relic densities are possible. 

Could the cosmological embarrassment be avoided by altering the composition of the 
LSP into which the gravitino is supposed to decay? The gravitino decay amplitude could be 
suppressed, and the lifetime estimate (56) correspondingly increased, if the LSP did not have 
photino component, for example if it were essentially a pure Higgsino. However, in view of 
the discussion above, in realistic models it seems difficult to suppress the (5 — y — 7 coupling 
enough to increase the gravitino lifetime sufficiently to relax the bound (57) adequately, 
bearing in mind the bounds applicable for longer lifetimes. 

Finally, we consider briefly the situation if R parity is violated. First we consider R- 
violating decays of the lightest MSSM particle, assumed to be the lightest neutralino y. As 
discussed earlier, its abundance is likely to be Y^ ^ 10 “^°, which conflicts with our bounds 
nnless ^ 10^ s Considering now the decays of the gravitino, assnming it to be heavier 
than X, there are two options to consider. The simplest possibility is that it decays in the 
same way as in the R-conserving case: (9 —y + 7 , etc., in which case the previous R- 
conserving gravitino analysis applies, and we again conclude that Tq ^ 10 ^ s. Alternatively, 
the dominant G decays might violate R parity. In this case, the bonnd Ty. ^ 10"^ s again 
applies, and, as discussed above, the same bound applies to the decays of the lightest MSSM 
sparticle. 

7 Discussion and Conclusions 

We have re-examined in this paper the upper limits on the possible abnndance of any nn- 
stable massive relic particle that are provided by the snccess of Big-Bang Nncleosynthesis 
calculations. A new aspect of this work has been the use of cosmic microwave background 
data to constrain independently the baryon-to-photon ratio, which was not possible in previ¬ 
ous studies of this problem. We have also incorporated in our analysis the an npdated suite 
of photonuclear and nuclear cross sections, and calculated both analytically and nnmerically 

®Note that the case ^ 10^ s requires further consideration, going beyond the scope of this paper. 
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the network of reactions induced by electromagnetic showers that create and destroy the 
light elements deuterium, ^He, ^He, ®Li and ^Li. 

It was pointed out in previous work that considerations of the deuterium abundance alone 
would allow certain exceptional regions of parameter space with relatively large abundances 
of unstable particles. However, as shown in this paper, considerations of the abundances of 
^He and ®Li exclude these particular regions. 

We have illustrated our results by applying them to massive gravitinos. If they weigh 
~ 100 GeV, their primordial abundance should have been ^ 5 x 10“^"^ x (100 GeV/m^), 
corresponding to a reheating temperature Tr < (1.9 — 7.5) x 10^ GeV. This could present a 
potential difficulty for some models of inflation and leptogenesis. We have discussed various 
scenarios for evading this potential embarrassment, for example by varying the gravitino 
mass, or by postulating an alternative scenario for baryogenesis, such as non-thermal lepto¬ 
genesis or the Affieck-Dine mechanism. 

This example of the gravitino illustrates the power and importance of the cosmological 
upper limits on the abundances of unstable massive particles. Extensions of this analysis are 
clearly desirable. For example, it would be valuable to combine our analysis of electromag¬ 
netic decay cascades with a similar analysis of hadronic showers, a topic that lies beyond the 
scope of this paper. 
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A Cross Sections 


Kawasaki, Kohri and Moroi [13] have provided a useful table of reactions and references to 
nuclear data. We have supplemented the data by using tabulations made by the National 
Nuclear Data Genter (NNDG) [ 68 ] and the NAGRE collaboration [69]. The relevant cross 
sections are listed here for convenience, and we refer the interested reader to the NNDG and 
NAGRE websites for further references on the data. 

We have computed thresholds and IQI values using the mass data of Audi and Wapstra[70], 
available on the US Nuclear Data Program website [71]. Reverse reaction data were some¬ 
times available, in which cases we used detailed balance to transform the data into forward 
data. The equations of detailed balance for the reactions 7 -I-T A + B and P-l-T ^ A + B 
are: 


0''1+T^A+B 

(^P+T^A+B 


QaQb 




(^A+B^T+^f 


(1 + ^AB)gT \ 

(1 + SpT)gAgBmAmBEcra(A, B) 

(1 + 5AB)gpgTmpmTEcm{P. T) 


(63) 

(64) 



where the gi are the statistical weights of each species, /i is the reduced mass of the system 
A + B, and Ecm{x,y) is the center-of-mass energy of the system x + y. See Blatt and 
Weisskopf [72] and Fowler, Caughlan, and Zimmerman [73] for discussions on these relations. 

In the numerical hts, all energies are in MeV. In a few cases, which we have noted, the 
hts are those previously published. Otherwise, we have adopted a specihc empirical form 
for the nonresonant parts of the cross sections. This form is the product of a power law in 
photon energy E^, and a power law in photon energy above threshold, E^ — |(5|. We have 
found that expressions of this type provide a simple but accurate representation of the data. 


1. d( 7 ,n)p E^^th = \Q\ = 2.224573 MeV [74]. 


<j{Ey) = 18.75mb 


E 




E 


2. t( 7 ,n)d E^^th = \Q\ = 6.257248 MeV [75, 76]. 




1.65 


E3.6 


3. t( 7 ,np)n 7^7,= |Q| = 8.481821 MeV [76]. 
\Q\^-%E, - |g|)2-3 


<j{E^) = 26.0mb- 




4. 3He(7,p)d Ey^th = IQI = 5.493485 MeV [77, 78]. We use reverse reaction data from 
NACRE [69, 79].’ 


(t{E^) = 8.88mb 


E3-4 


1.65 


5. 3He(7,np)p E^^th = IQI = 7.718058 MeV [78, 80]. 


<j{Ey) = 16.7mb 


\Qr^{E,-\Q\f 

^4.25 


6. ^He( 7 ,p)t E^^th = |g| = 19.813852 MeV [81, 82, 


^{Ey) = 19.5mb 


iQ^Ey-m 


1.0 
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7. ^He(7,n)3He E^^th = \Q\ = 20.577615 MeV [82, 84, 85, 86, 87], 


(t{E^) = 17.1mb 


W-\E,-\Q\Y-^ 


^He( 7 ,d)d E^th = \Q\ = 23.846527 MeV 
from NACRE [69,’90]. 


We use reverse reaction data 


(y{EY) = 10.7mb 


\Qr\E,-\ Q\) 


3.4 


9. ^He( 7 ,np)d E^^th = \Q\ = 26.0711 MeV [84]. 

\Qr{E, - |g|)3-° 


<j{EY = 21.7mb- 


E7-0 


10. ®Li( 7 ,np)^He E^^th = \Q\ = 3.698892 MeV [91, 92]. 


<j{EY = 104.mb 


\QnE,-\ Q\) 

E7-0 


4.7 


11. ®Li(7,X)3A E^^th = \Q\ = 15.794685 MeV [93]. 


3.7 exp 


1 rE^ - 19.0' 
3.5 


+ 2.75 exp 


1 rE^ - 30.0 
3.0 


+ 2.2 exp 


1 rE^ - 43.0 
5.0 


12. ^Li( 7 ,t)^He E^^th = IQI = 2.467032 MeV Ecm = E^ — |g| [92]. We use 
reverse reaction data from NACRE [69], with modifications from Cyburt, Fields, and 
Olive [16, 94]. 

nin. ,.^2371MeV2\ f-2.595A\ 

<j{EY = 0.105mb I- — -) exp I —I x 


V J 


VK 


exp 


(-2.056Ee™) (l. + 2.2875^2^ - 1.1798^3^ + 2.5279E,^J 
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13. ^Li( 7 ,n) 6 Li E^^th = |Q| = 7.249962 MeV [95, 96, 97]. 


. 0.176mbM7(47K^ + 1205mb + 


0.06mb 




E10.5 


1 _|_ / i?em—7.46 y 

V 0.188 J 


14. ’^Li(7,2np)"^He E^^th = |Q| = 10.948850 MeV [92]. 


(y{E^) = 122.mb 


\Qr{E,-\Q\) 

E]-^ 


3.0 


15. ’^Be( 7 ,^He)^He E^^th = |Q| = 1.586627 MeV Ecm = E^ — |(5|. We use 
reverse reaction data from NACRE [69], with modifications from Cyburt, Fields, and 
Olive [16, 98]. 


uiE^) = 0.504mb 


/2371MeV2\ /-5.1909\ 


exp 


'\/ Ecm ) 


X 


V J 

exp (-0.548Ee„7) (l. - 0.428^^^ + 0.534^3,^ - 0.115E; 


4 

cm 


16. ^Be( 7 ,p)®Li E^^th = |Q| = 5.605794 MeV. We use reverse reaction data from 
NACRE [69, 99]. 


<j{Ey) = 32.6mb 




^ 12.0 


+ 2.27xl0‘’mb^ 


Z7'21.8335 

7 


17. ’^Be(7,2pn)"‘He = \Q\ = 9.30468 MeV. No data exist. We assume isospin 

symmetry and use the data for the reaction ^Li( 7 , 2np)^He [92]. 


<j{E^) = 133.mb 


\Qr{E,-\Q\) 

E7-0 


3.0 


4rT„/'4- 

reverse reaction data [100]. 


18. "^He(t,n)®Li Ep^th = 8.386972 MeV \Q\ = 4.78293 MeV Ecm = mT+mt Ep- 


. X A{Ecm-\Q\)/\Q\Y^ A{Ecm-\Q\)/\Q\A^ A{Ecm-\Q\)/\Q\A 

a{E^) = 1940.mb^V^^h^—+ 46.2mb ^^ , ^ 32.8mb-^^ 


^ _|_ |' gcm-5.03043 y 




1 I ( Ecm-7.5239 y 
^ \ 25.0 J 


19. 4He(3He,p)6Li Ep,th = 7.047667 MeV \Q\ = 4.019167 MeV Ecm = We 

use reverse reaction data [ 101 ]. 

ciE-,) = lT0.mbI<^^MMa7 ^ + 32.1mbF“‘ " 


^ I' Ecm-5.8192 y 


^ |' Ecm-7.5192 y 


^ |' Ecm-12.019 y 


31 




Center—of—Mass Energy (MeV) 


Figure 9: Cross-section data for the reaction ‘^}le(t,n)^Li are plotted versus the center-of- 
mass energy. The solid curve is the parameterization we use, whilst the dashed curve is that 
adopted in previous studies [12, 13]. 
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^He(a,p)®Li 



Center—of—Mass Energy (MeV) 


Figure 10: Cross-section data for the reaction ^He(^He,]9/Li are plotted versus the center- 
of-mass energy. The solid curve is the parametrization we use, whilst the dashed curve is 
that adopted in previous studies [12, 13], 
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